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Abstract. Let A = Ei x E2 be be the product of two elliptic curves over Q, 
both having a rational five torsion point Pi. Set B = A/ {{Pi, P2)) ■ In this 
paper we give an algorithm to decide whether the Tate-Shafarevich group of 
the abelian surface B has square order or order five times a square, assuming 
that we can find a basis for the Mordell- Weil groups of both Ei , and that the 
Tate-Shafarevich groups of the Ei are finite. 

We considered all pairs {Ei, E2), with prescribed bounds on the conductor 
and the coefficients on a minimal Weierstrass equation. In total we consid- 
ered around 20.0 million of abelian surfaces of which 49.16% have a Tate- 
Shafarevich group of non-square order. 



1. Introduction 

Let A be an abelian variety over a number field K. Then the Tate-Shafarevich 
group III{A/K) plays an important role in understanding the arithmetic of A. 
For example, it contains information on the tightness of the upper bound on the 
Mordell- Weil rank obtained by m-descent. Moreover the order of this group, which 
is conjectured to be finite, plays a role in the Birch and Swinnerton-Dyer conjecture. 

The Tate-Shafarevich group comes with a pairing, the so-called Cassels-Tate 
pairing, which depends on the choice of a polarization X : A ^ A"^: 

{; •)> : m{A/K) X m{A/K) Q/Z. 

Let II1{A/ K)nd denote the Tate-Shafarevich group modulo its maximal divisible 
subgroup. If A is an isomorphism, i.e., A is principally polarized, then the induced 
pairing on in(A/i(r)nd is non-degenerate. If moreover this pairing is alternating, 
then for all primes p the cardinality of the p-divisible part U1{A/ K)a(i[p°°] is a 
perfect square, thus if Ill{A/K) is finite, its order is a perfect square. 

Tate [17] showed that if A is an isomorphism and is also induced from a K- 
rational divisor on A then the Cassels-Tate pairing is actually alternating, as for 
example for elliptic curves. However, if dim A > 1 then A may not admit a principal 
polarization and even when A is principally polarized then this polarization need 
not to be induced by a i^T-rational divisor on A. Poonen and Stoll [TU] in fact showed 
that there exist genus 2 curves C such that #III( J(C)) is twice a square. Moreover, 
they showed that if one assumes that III(J(C)) is finite for all genus 2 curves C/Q 
then the density of the Jacobians of genus 2 curves that have non-square order 
Tate-Shafarevich groups exists, and they showed that it is about 13%. 
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For arbitrary abelian varieties Flach showed that if ^m{A/K) = kn^, with 
k square free, then k divides 2 times the degree of any polarization on A. Hence for 
principally polarized abelian varieties one has that #111 (^/X) is either a square 
or twice a square, if it is finite, but for general abelian varieties there are more 
possibilities. Stein |15| constructed for every prime number p < 25, 000 an example 
of a p — 1-dimensional abelian variety Ap/Q such that ^UI{Ap) = pn? . 

We restrict now to the case of dim A = 2. The constructions of Poonen-Stoll and 
of Stein yield examples of abelian surfaces such that ^'U1{A/K) is a square, twice a 
square or three times a square. One might wonder which further possibilities occur. 
Recently, the first author [4J showed that there exist abelian surfaces such that the 
Tate-Shafarevich group has order five times a square and seven times a square. 

In this paper we will have a closer look on the construction of abelian surfaces 
with Tate-Shafarevich group of order five times a square. The examples of (4] 
are members of a two dimensional family of abelian surfaces with a polarization 
of degree 5^. Moreover, one can show that for a general member of this family 
every polarization it possesses has degree a multiple of 5, thus they are not a priori 
excluded by Flach's theorem and might have a Tate-Shafarevich group of order five 
times a square. 

The construction of this family goes as follows. Let {E, O) be an elliptic curve 
over Q with a point P of order 5, then there exists a c? G Q* such that ((i?, O), P) 
is isomorphic to {{Ed, O), (0, 0)) with 

Ed ■ y + {d + l)xy + dy = + dx^ . 

Take now two numbers c?i, ^2 G Q* and consider Bd^^d^ '■— Ed^ xi^^s/ ((0, 0) x (0, 0)). 
Then A^i.da ■= Ed^ x Ed2 ^di.da is an isogeny of degree 5. Moreover, if the 
two elliptic curves are not isogenous, then all polarization on Bd^^d^ have degree 
divisible by 5. The i3di,d2's are the family we consider. In our case we know that 
III(^di,d2/Q) has square order, if it is finite, since it is isomorphic to the product 
of the two Tate-Shafarevich groups of Ed^ and Ed^- 

The behavior of the Tate-Shafarevich group under isogenics is well-known. This 
behavior is part of Tate's proof of the invariance of the Birch and Swinnerton-Dyer 
conjecture; for more on this see Section [2] The upshot of this is the following: 
Let (f : A ^ B he dju isogeny and assume that either ^'U1{A/K) or ^'U1{B/K) 
is finite (which implies that both are finite). Denote by (^^ : B^ — ^ A"^ the dual 
isogeny. For a field L D K denote by (pL ■ A{L) B{L) the induced map on 
L-rational points. Let S' be a finite set of places containing the primes where A has 
bad reduction, the infinite places and the primes dividing the degree of tp. Then 
the following holds: 

#III(A/iir) #ker coker(^^ -p-r #coker(/?/f^ 
#m(B/i^) " #ker<#coker(^^f j-^l #ker(p;f„ 

In Sections S] and [S] we show that for our choice of abelian surfaces the above 
mentioned cardinalities of kernels and co-kernels can be determined, provided one 
has a basis for the Mordell-Weil group of both Ed^ and Ed^ ■ (Actually something 
weaker is enough, see end of Sectional) Hence, given a basis for the Mordell-Weil 
groups of both elliptic curves we can determine whether #IH(_B/Q) is a square or 
a non-square. 

For all pairs {di,d2) with di — Ui/vi, such that max(|wi|, is bounded by 
N = 50,000 and the conductor of Ed^ is bounded by C = 10^, we computed this 
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product of cardinalities of kernels and cokernels. There are 2,445,366 such pairs 
and 47.01% of these surfaces have a Tate-Shafarevich group of non-square order. 
Also we computed these cardinalities for all pairs (di,c?2), such that the absolute 
value of the numerator and denominator of di is bounded by A'^ = 100. There 
are 18, 522, 741 of such pairs and 49.31% of them have a Tate-Shafarevich group of 
non-square order. Based on our computations we expect that the density of abelian 
surfaces Bd^^d2 with non-square Tate-Shafarevich groups exists and is around 50%. 
For some heuristics see the end of the last section. 

The outline of this paper is as follows. In Section [2] we discuss some preliminar- 
ies and in Section [3] we explain in more detail the construction of the considered 
familiy of abelian surfaces. In Section |4] we discuss how we can calculate the global 
quotient and which conditions on Ed-^ and Ed2 are needed for this. In Section [5] 
we discuss how we calculate the local quotient, which turns out to be much easier 
computationally. In Section IH] we sketch the algorithm used for the computations 
of the densities and finally in Section [7] we discuss the obtained results. 



Let K he a. number field, and let Gk be the absolute Galois group Gal{K /K). 
For a (finite or infinite) place v oi K denote by K^, its completion with respect to 
V and Gk^ its absolute Galois group. 

Let A/K be an abelian variety. Denote by the dual abelian variety. Then 
the Tate-Shafarevich group of A/K is defined as 



where the product is taken over all finite and infinite places of K. Let ip : A ^ B 
be an isogeny of abelian varieties, then the (/j-Selmer group of A/K is defined as 



Here [•] means "kernel of". 

The Tate-Shafarevich group is torsion. It is conjectured to be finite and the (f- 
Selmer group is known to be finite. The m-torsion subgroup of the Tate-Shafarevich 
group fits in an exact sequence 



I.e., it measures the difference between the m-Selmer group and A{K)/mA{K). 
In theory the m-Selmer group is computable, hence the Tate-Shafarevich group 
measures the difference between the upper bound on the Mordell-Weil rank obtained 
by doing m-descent and the actual Mordell-Weil rank of A. 

The Tate-Shafarevich group plays also a role in the Birch and Swinnerton-Dyer 
conjecture: 

Conjecture 2.1 (Birch and Swinnteron-Dyer). Let A/K he an abelian variety and 
L(A, s) its L-series. Set r := T]<iA{K). Then Iil{A/K) is finite, L{A, s) has a zero 
of exact order r at s — I, and 



2. Preliminaries 





^ A{K)/mA{K) ^ S^'"^{A/K) UI{A/K)[m] 0. 



L{A,s) _ T#Ul{A/K)RAUlAiK.^H 

™1 (s - 1)'- " #A{K)tor#A''{K)tor 
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The left hand side of this conjecture is invariant under isogeny. Cassels fl] 
(dim A = 1) and Tate [ITj (diniv4 > 1) proved that the right hand side is also 
invariant under isogeny. I.e., if 1^9 : ^ _B is an isogeny then 

#m{A/K) _ i?B#A(j^)to.#A^(jr)torn/B(gJ l^k 
#m(B/i^) " i?^#i3(i^)tor#S^(i^)tornL(K„)l^l-' 

This formula was used by Schaefer and the second author to provide examples of 
elliptic curves with large Selmer groups, by Matsuno [S] and by the second author 
ini to provide examples of elliptic curves with large Tate-Shafarevich groups and by 
Flynn and Grattoni |3J to compute several Selmer groups. 

However, for calculation purposes the right hand side is not suitable. One can 
rewrite the right hand side as follows: For a field L Z) K lei ipi^ denote the group 
homomorphism ifL ■ ^{L) B{L). Then 

^UI{A/K) #kcr(^/f#coker(^^ -p-i- ^cokenpK^ 
#m{B/K) = #kcr(^^,#coker(pK 11 #ker(^K^ ' 

We will call the first factor with the ipK the global factor, the second factor with 
the ifK^ the local factor. If is a finite prime of good reduction and v does not 
divides the degree of the isogeny then # coker (pK„ = # ker ipK^ , hence the product 
on the right hand side is a finite product, where only the bad primes, the infinite 
primes and the primes dividing the degree of the isogeny are taken into account. 

It is known that if the analytic rank of an elliptic curve is at most 1, then its 
Tate-Shafarevich group is finite and the analytic rank equals the Mordell-Weil rank; 
otherwise we will assume these two conjectures. 

3. Constructing a family of abelian surfaces 

We will construct a two-dimensional family of abelian surfaces B/K , whose mem- 
bers are quotients of products of two elliptic curves Ei , E2 by an isogeny of degree 
5. Therefore #U1{B/K) ■ 5° = #lll{Ei x E2), for an a G Z. Since #IU{Ei x E2) 
is a square it follows that ^111(5//^) modulo squares is one of {1, 5}. Additionally 
we have that for a general member of this family every polarization has degree 
divisible by 5. Thus Flach's theorem does not restrict us further. 

Let G/K be a group scheme of prime order £. Let Ei,E2/K be two elliptic 
curves such that G is a subgroup scheme of both Ei and E2 ■ Let A = Ei x E2 and 
B = A/G. Then ip : A —i' B has degree £. Moreover, one can show that either Ei 
and E2 are isogenous or every polarization on B has degree a multiple of £. Hence 
for general Ei, E2 we are in the second case. 

Consider the case G — Z/£Z, i.e., G is generated by a if -rational point. Since 
for ^ > 4 the functor Yi{£) is representable one has a universal family of elliptic 
curves E with a point P of order £. In the case £ — 5 the universal family is given 

by 

Ed:y'^ + {d+ l)xy + dy = + c^a;^ P = (0, 0), 
for any d G K* with d^ + lid — 1 7^ 0. The four non-trivial 5-torsion points are 
(0, 0), (-d, d^), (-d, 0), (0, -d). If we move (0, -d) to (0, 0) and bring the curve in 
standard form we obtain Ed- If we move (— d, d^) or (— d, 0) to (0, 0) and bring the 
elliptic curve in standard form then we obtain E^i/^- 

We restrict now to the case Tf = Q, ^ = 5, and G is generated by a Q-rational 
point. Fix di and d2 in Q* , set A := Ed^ x Ed^ ■ The rational 5-torsion subgroup of A 
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has four diagonally embedded subgroups of order 5. Let G = Z/5Z be one of those, 
i.e., it is the subscheme of A generated by (0,0) x [n](0,0), with n £ {1,2,3,4}. 
Let B := A/G. Then S is a candidate for an abelian surface such that in(i3/Q) 
has order five times a square. To actually check whether Ill(i?/Q) has non-square 
order we will now calculate both the local and the global factor. 

Note, that the 16 surfaces B/Q one obtains by replacing di by —l/di and using 
the four values of n break into two pairs of 8 isomorphic surfaces. For fixed di,c?2 
the surfaces corresponding to n = 1,4 lie in one of these isomorphism classes and 
those for n = 2,3 in the other one. We will see in the next two sections that for 
fixed di,d2 the size of III(S/Q) is independent of n, thus all 16 surfaces will have 
Tate-Shafarevich groups of same cardinality. Therefore, for the computations we 
will only consider the case rfi , ^2 > and n — 1. 

Let A' be the quotient of Ed^ x Ed^ by ((0,0) x 0,0 x (0,0)) and E'^^ be the 
quotient of Ed^ by ((0,0)). The isogeny A A' factors as A^ B ^ A'. Consider 
now the dual picture 

{A'y -^B^ ^ A"". 

Since A and A' are products of elliptic curves, they are principally polarized. There- 
fore we have the following factorization 

A' ^B"" ^ A. 

The kernel oi A' A is Cartier dual to the kernel of A — >■ A', and hence is 
isomorphic to (ms)^. The kernel of A' — >• is isomorphic to /j.^ embedded with 

(1, -n) in (/^5)2. 

Summarizing we have the following diagram: 




Lemma 3.1. Suppose L = Q. Then kcr^sq = Z/5Z and ker(/9q = 0. 

Proof. Since A[ip] = Z/5Z it follows that ^'[<^^] = /iXg. Taking Q-rational points 
yields the lemma. □ 

Lemma 3.2. Suppose L = R. Then kenpR = Z/5Z and cokeri^R, = 0. 

Proof. The first assertion is automatic. The non-trivial element in Gal(C/R) acts 
on the fiber of an clement of i3(R) under ipc either by swaping elements or fixing 
them. Since the degree of (f is not divisible by 2 at least one element in the fiber is 
fixed, hence lies in A(R). □ 

Let S be the set of primes where A has bad reduction, together with 5. Using 
the above lemmas it follows that 

#m{A/Q) ^ #cokery)^ #cokery.Q^ 
#m(i?/Q) " #cokcr^Q i^i #ker^Q„ ' 

In the next two sections we will first explain how to determine the first factor, the 
global factor, then how to determine the second factor, the local factor. 
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4. Determining the global factor 

To determine ^ coker assume for the moment that one has a basis for the 

Mordell-Weil groups (Q), -Bd^ (Q), ^^rf^ (Q) and E'^^{q,). We will now explain 
how one can determine coker (/?q and coker (/?q from this. 

Using = ip^ o^^ we obtain a surjective homomorphism coker — cokeriy9q. 
With Hilbert's Theorem 90 we obtain 

The surjection coker pq — )■ coker iy9q becomes {x,y) a;"/y as the map from 
(Q*/Q*^)^ to Q*/Q*^. One sees immediately that the image of this map is inde- 
pendent of n, hence to compute coker (/?q we may set n = I. In order to determine 
coker iy9q it suffices to determine a basis for both coker 77^'q in Q*/Q*^. This can 
be done quite easily following [14] Exercise 10.1]: Suppose that / is a function on 
with divisor 5(0,0) — 50. Then there exists a unique constant c G Q*/Q*^ 
such that the map 

coker^q ^ Q7Q*^ 

sending P ^ (0,0), O to c/(P) mod Q*^, is a well-defined and injective group 
homomorphism and its image agrees with the image of the natural embedding 
of coker r^j^q into H^{Gq, E'^, [??/]) — Q*/Q*'^. In our case we can take the function 
/ = —x^ + y + xy and the constant c = 1. The point (0, 0) is mapped to d^^ and 
O to 1 by linearity. 

An element of Q*/Q*^ is determined by the valuations at each prime. Write 
now d = u/v and let S be the set of all primes p dividing five times the minimal 
discriminant of Ed, i.e., p \ buv{v? + Wuv — v^). Define 

Q(S', 5) := {x e Q7Q*^ | Vp{x) = mod 5, Vp i S). 

From the same exercise from [T3] it follows that / (coker yyq) C Q(S', 5). Hence we 
can represent an element of coker r^q by its valuation at each prime number p ^ S. 
Once the cokernels of both tj^q are established, the cokernel of can be computed 
easily. 

To determine the cokernel of ipQ we use the following exact sequence 

— ?> ker(?/;Q)/<^(kerpQ) — > coker (^q coker pq coker i/jq -> 0. 

Note that kei-{ipQ) = (^(kerpq). Set K := Q(C5), where ^5 is a primitive fifth root 
of unity. Then the restriction map 7?^(Gq,Z/5Z) {G k , 'Zi / 5Z) is injective 

since the kernel of this map has exponent dividing both : Q] = 4 and #Z/5Z. 
Since A[ip], A[p\ and Blip] are isomorphic to fi^, fi^ x /Xg and fi^ over K we obtain 
the following commutative diagram with embeddings as vertical maps. 

*- coker (pQ *- coker pq *- coker ipQ 



*- K*/K*^ *- {K*/K*^)^ K*/K*^ *- 
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As above the lower second horizontal map is just (x, y) i— t- a;"/?/. Hence to determine 
cokeriy9Q it suffices to determine the kernel of x"'/y on coker77i^Q x coker7]2,Q 
coker^Q. Again this is independent of n. We do this as follows: 

(1) For some d Q K there is a /^-isomorphism t : E'^^ E^^, sending a generator 
of kerr;^ to (0,0). The map f : E'^ K*/K*^ is then P ^ -x{T{P)f + 
y{T{P)) + x{T{P))y{T{P)). Hence we have to determine t. This can be 
done easily for each individual curve 

(2) To represent elements in cokerr/q C K*/K*^ note that the class number 
oi K* equals 1. Set 

K{S, 5) {x e K*/K*^ \ v^{x) = mod 5, Vp ^ S], 

where S contains all primes p of X being a bad prime of Ed or dividing 5, 
i.e., all primes p of X lying over a primes p, such that p \ 5uu(m'^ + 11uw — u^). 
From [m Exercise 10.9] it follows that /(cokerr^q) C K{S,'o). Hence to 
represent elements in coker Tyq we have to fix a generator tp for each prime 
p G 5, and we have to fix generators for the unit group of K modulo fifth 
powers. The field K is well-known and it is easy to see that the unit group 
is generated by —C,^ and (1 Q^). Hence we can write 

/(p)^c5""(i + C5r nc^^^''^^ 

pes 

modulo fifth powers. 

Remark 4.1. We can weaken the assumption of having a basis for the Mordell- 
Weil groups ^E^^ (Q), i?^^ (Q), i?^ (Q) and E'^ (Q). It is actually sufficient to just 
have generators of a finite index sublattice of these four groups, such that the index 
is not divisible by 5, i.e., the generators of infinite order are not divisible by 5 
modulo torsion. This is the case, since the images of such sublattices in the co- 
kernels of ?7j^, respectively rji, are the complete co-kernels. Also it is sufficient to 
just know such sublattices of i?dj(Q) and £'d2(Q), since suitable dual sublattices 
can be easily computed using the isogenies rji. One only has to calculate the images 
of the generators under rji and then check if their span contains points divisible by 
5 modulo torsion. 



5. Determining the local factor 

We want to calculate '^^"^^-cip^^ ^^-^ ^^"-^ primes p and for p = deg = 5. Since 
the kernel of i^q is generated by a Q-rational point it follows that ^ker^sq = 5. 
The size of the co-kernel of t/jq^^ depends on the reduction of Ed^ and Ed2 , but turns 
out to be independent of n. 

For rj :— rji, we first describe how coker Tyq^ depends on the reduction type of 
E :— Edi- Write di —: u/v with 7i,w G Z and gcd(u,w) = 1. Then E has the 
following global minimal equation 

E : y'^ + {u + v)xy + uvy = x^ + uv^x"^ 

and discriminant — (ww)^(w^ -|- lluv — v^). 

Lemma 5.1. E has the following reduction type at a prime p: 

(1) If p \ uv then the reduction is split multiplicative and the point (0,0) does 
not lie on the identity component of the Neron model of E. 



8 



STEFAN KEIL AND REMKE KLOOSTERMAN 



(2) If p I + lluv — then (0, 0) lies on the identity component of the Neron 
model of E and either p — 5, or p = ±1 mod 5 holds. If p = 5 the reduction 
is additive, if p = 1 mod 5 then the reduction is split multiplicative, and if 
p = 4 mod 5 then the reduction type is non-split multiplicative. 

Proof. Let E he E mod p and E^s be the smooth locus of E. If p | uv then E has 
equation + axy — x^, for some non-zero a G Z/pZ. In particular, (0,0) mod p 
is a node of E and the tangent cone is generated by a; = —ay and y = 0, hence the 
reduction is split multiplicative. Since (0, 0) reduces to the singular point of E this 
point does not lie on the identity component of the Neron model of E. 

If p I + lluv — v"^ then the reduction of (0, 0) is both on i?ns and is non-trivial. 
In particular the order of the reduction of (0,0), which is 5, divides #i?ns(Fp). If 
the reduction is split multiplicative then this group has order p—1, if the reduction 
is non-split then this group has order p -I- 1, and if the reduction is additive then 
this group has order p, i.e., p = 1 mod 5, p = —1 mod 5, and p = 5 respectively. □ 

Let E' := E'^, be the isogenous elliptic curve. Denote by ce,p and ce',p the local 
Tamagawa numbers, i.e., the number of components of the Neron model. 

Lemma 5.2. For the Tamagawa quotient we have 



Proof. Since rj has degree 5 it follows that that — 5" for some a G Z. If the 
reduction is different from split multiplicative then c^; ^ and CE'.p are at most 4, 
hence a — and c^ ^ = CE'.p. 

In m Proposition 2.16] it is shown by using Tate curves that if the reduction 
is split multiplicative then a G { — 1, 1}, depending on whether the kernel is on the 
identity component of the Neron model or not. □ 

If p f degT? = 5 then from [Til Lemma 3.8] it follows that *^tl7nQ'' = 
Using this it follows easily that 

Lemma 5.3. Suppose p is a prime different from 5. 

(1) If p I + lluv — v"^ and p = 4 mod 5, then cokerr/Q^ = Z/5Z. 

(2) If p I -I- lluv — v'^ and p = 1 mod 5, then cokerr^Q^ = (Z/5Z)'^. 

(3) If p I uv, then cokerr/q^ = 0. 

(4) If p is good for E, then coker?7Qp — 'L/^'Z. 

Now coker?7Qp C II^{Gq^,Z/5Z). From [H Section II. 5 Theorem 2, Proposi- 
tion 17] it follows that for p \ deg 77 = 5 



with a = 1, if p = 4 mod 5, and a = 2, if p 1 mod 5. From this it follows that 

Proposition 5.4. Suppose p is a prime of bad reduction for E and p 7^ 5 with 
p \ u'^ + lluv - v'^ . Then cokerT^q^ = H^{Gq^,Z/5Z). 

We now return to our abelian surface A. Then the above proposition enables us 
to determine cokeri/jq^ for bad primes different from 5. 




5, if p I -|- lluv — v'^ and p — 1 mod 5, 
1, otherwise. 



#H\Gq^,Z/5Z) = #i/"(Gq^,Z/5Z)#iJ"(Gq^,/x5) = 5' 



a 
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Proposition 5.5. Suppose p is a prime of bad reduction for A, p ^ 5. Then 
cokeriy9Q^ is isomorphic as an abelian group to 

0, if P \ U1V1U2V2, 

(Z/5Z)^, if p I gcd(M^ + lluiwi — vl,ul + llti2W2 - ), P = 1 mod 5, 
Z/5Z, otherwise. 

Proof. Recall that coker(pQ^ — ker (cokerr^i.Qp x coker772.Qp — >■ coker'i/^Qp) , which 
equals 

(cokerr/i,Q^ x cokerr/a.qj n ker (i/i(GQ^, Z/5Z)2 ^ H^Gq^, Z/5Z)) . 

The surjective map H^{Gq^,Z/5Z)'^ -> H^{Gq^,Z/5Z) is given by (x, y) i-> nx—y. 
Suppose that p | U1V1U2V2, then by Lemma l5.3l we have that coker?7i.Qp = 0, for at 
least one i, and therefore cokeiipQ^ = 0. 

Suppose now p -j" uiViU2V2- By assumption one of the E^^ has bad reduction 
at p, let's say Ed-^. Since p \ 5uiVi it follows from the above proposition that 
cokeryyi^Qp = H'^{Gq^,Z/5Z) and hence cokeiipQ^ = coker 772,Qp. Now E^^ has 
either additive or good reduction. The reduction of Ed2 is additive if and only if 
p I gcd(wf + llwiWi — vl,U2 + IIM2W2 ~ ^2). Now apply Lemma [5.31 to deduce the 
structure of coker?72,qp, hence the structure of cokeriySQ^. □ 

It remains to check the case p = 5. As before, we first have a look at the elliptic 
curve E. If 5 | uv then as above the reduction is split multiplicative and = ^■ 
Using Tate curves one easily shows that coker?7Qp = 0. 

If 5 I + lluv — then the reduction is additive. In particular the component 
groups of E and E' have the same order, which is also the case if the reduction 
is good. Therefore — 1- The isogeny rj : E ^ E' can be written as a power 
series in one variable in a neighbourhood of the point O. Again from [111 Lemma 
3.8] it follows that 

#ker?7Q5 

where |?7'(0)|5 is the normalized 5-adic absolute value of the leading coefficient of 
the power series representation of 77 evaluated at 0. This can be easily computed 
using Vein's algorithm [TB]. In [H Lemma 4.1, Proposition 4.2] it is shown that in 
the additive case vsiu"^ + lluv — v^) G {2, 3} and that if v^^u? + lluv — w^) — 2 
then |7?'(0)|5 = 1, and if v^iu^ + lluv - w^) = 3 then |7?'(0)|5 = 1/5. If E has good 
reduction at p = 5 then it follows that # coker rjQ — ^ ker 77Q , because in this 
case we also have |77'(0)|5 = 1. We summarize as follows. 

Lemma 5.6. Suppose p — b. 

(1) If p I uv, then cokerT^q^ = 0. 

(2) If p is good for E, then coker ?7Qp = Z/5Z. 

'(Z/5Z)2, 53 I u2 + lluv - 
Z/5Z, otherwise. 

Now we can calculate coker (pq^ in the remaining case p = 5. 

Lemma 5.7. The cardinality 0/ coker (pq^ equals the cardinality ofkenpQ^, unless 

(1) 5 I U1V1U2V2. In this case # coker i^q^ — 0. 

(2) 5'^ I gcd(?i^ + lluiUi — wj, U2 + IIU2W2 — ^^2)- In this case ^ cokei ifiQ^ = 5'^ . 



(3) If p \ u^ -\- lluv — v"^ , then coker jyq^ = 
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Proof. If cokcry^i^Qj; = 0, for one i, then cokenpQ^ — 0. The first condition is 
equivalent with 5 | uiViU2V2- 

Suppose now that cokeryyi Q^ ^ for both i, which imphes that p = 5 is ad- 
ditive or good for Ed^- We need two facts from [T31 Section II. 5 Proposition 18, 
Theorem 5], namely ^H^{Gq^,Z/5Z) = (Z/SZ)^ and H^,^{Gq,,Z/5Z) = Z/5Z. 
As in the previous proposition we have that if cokerr/i^q^ = H-^{Gq^,Z/5Z), 
then cokeriySQg = cokerr/2,Q5 and vice versa. This gives the second case of the 
lemma, since coker77i.Q5 ~ (Z/5Z)^ if and only if 5^ | w| + lluiVi — vf, and 
coker?7i_Q5 = Z/5Z otherwise. 

It remains to consider cokerT^i Q^ — coker772.Q5 = (Z/5Z). In this case one 
can show that cokerT^i.Qj = iJ,j^(GQ5 , Z/5Z), for both i; see [U Proposition 2.10, 
Proposition 3.5] and [12, Section 3]. Thus the kernel of cokerr^i.Qg x coker?72,Q5 — > 
cokerV'Qj, which equals cokercpq^, has five elements. This finishes the proof. □ 

Putting everything together yields 

Proposition 5.8. Let p be a prime. Then 

# cokei ipQ^ 
# ker ipQ^ 

is a non-square if and only if one of the following occurs 

(1) p I U1V1U2V2, 

(2) p I gcd(wf + lluiwi — vf, U2 + IIU2W2 — w|) and p =1 mod 5, 

(3) p^ I gcd(M2 + lluiwi — v1,U2 + llM2f2 " v|) and p — 5. 

6. Algorithm 

The code was implemented in Sage [16] and is available at [5j. The algorithm 
consists of two main steps and an initialization step, which we call step 0. In step 
1 one creates a database of elliptic curves having a point P of order 5, which are 
parametrized by two coprime positive integers (u, v). One has to specify which pairs 
(m, v) one wants to consider. In step 2 one takes such a database of elliptic curves 
Ed, for d — u/v, and goes over all pairs of these curves and determines whether the 
order of the Tate-Shafarevich group of the abelian surfaces B = Ed^ xEd2/ {{Pi, P2)) 
is a square. For trivial reasons, pairs of the same elliptic curve are omitted and 
pairs are considered to be without order. 

Step 0: Fix a (large) integer M . For each prime number p < M determine the 
prime ideals p oi K = Q(C5) above p and fix an ordering of them. Then fix for each 
prime ideal p a generator tp . 

Step 1: Fix a positive integer N. For each pair of coprime positive integers (u, v), 
such that max(u,t;) < N, we collect the following data associated to E Ed, for 
d = u/v. (Optional: Filter the pairs (u,f) by other limitations, e.g., considering 
only those for which the corresponding elliptic curves have conductor < C.) 

• Collect all the primes dividing 5uv{u^ + lluv — w^) in a set S. 

• Collect all the primes dividing uv in a set T. 

• Collect all the primes p = 1 mod 5 dividing + lluv — v'^ in a set U. 

• If W5(u^ + lluv — v'^) — 3, put also p — 5 into the set U. 

• Determine the analytic rank r of E. 

• Determine a system of r generators of a sublattice A of E{Q), such that 
the points of infinite order modulo torsion are not divisible by 5. Take the 
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image of A in Q{S,5) to determine a basis P of cokerr^g C Q{S,5). The 
data for each basis element consists of a pair for each prime in S, where 
the first entry is the corresponding element in S and the second entry is 
the exponent as an element in Z/5Z. 

• Calculate the image of A under rj in -E"(Q) and determine which image 
points are divisible by 5 modulo torsion. Divide if possible and determine 
the non-trivial 5-torsion points of E'{Q,) to get a sublattice A' of £"(Q), 
such that the points of infinite order modulo torsion are not divisible by 5. 
Use this information to get dim coker rjQ . 

• Take the image of A' in K{S,^) to determine a basis Q for coker jjq c 
K{S, 5). The data for each basis element consists of a pair for each prime 
in S and a pair for the units, where the first entry is the corresponding 
element in S, respectively 1, and the second entry is a list of elements 
in Z/5Z, which contains as many entries as there are prime ideals p in 
K over p, respectively two entries, and these entries are the exponents 
corresponding to the prime ideals (tp) with the chosen order, respectively 
the exponents of the units. 

Step 2: To determine whether the order of Tll{Bdi,d2/Q) is ^ square, for each 
pair of pairs {ui,vi), (w2, V2) from the first step (modulo ordering and equality), do 
the following: 

. Set L := -#(Ti U T2) + #{Ui U2). 

• Fix an ordering for <S := U S'2. 

• Write out the elements from Pi U P2 into a matrix with respect to <S. This 
gives a matrix with entries in Z/5Z. Calculate the rank of this matrix, 
which equals the dimension of coker i^g . 

• Write out the elements from Qi U Q2 into a matrix with respect to the 
prime ideals (tp) lying over the primes of S (and with respect to the units). 
This gives a matrix with entries in T^/^T^. Calculate the rank of this matrix, 
which equals the dimension of coker i/'q • 

• Set G := dim coker (/Jq — dim coker t/i^q — dim coker 772, Q + dim coker V^q . 
(Recall that dim coker 9?q = dim coker ?7i^q + dim coker 7^2. q — coker t/jq.) 

Then the local factor (without the infinite prime) is a non-square if and only if 
L is odd, and the global factor (without the kernels) is a non-square if and only if 
G is odd. Since the contribution of the infinite prime and the kernels cancel, we 
have that III(i?di,d2/Q) non-square order if and only if L -|- G is odd. 

The constructed databases and obtained results are given in the next section. To 
conclude this section, we make some comments on the implementation. Step in 
the cases considered is not computational demanding. For example, on a desktop 
computer it may take some seconds up to a few minutes to compute all generators 
for all prime ideals lying over all primes up to 500, 000. Step 2 is also no problem. It 
consists only of simple set operations and calculating ranks of small matrices with 
coefficients in Z/5Z. Even a few million of pairs of elliptic curves can be considered 
in under an hour. 

The computational demanding part is step 1. There are two main issues. The 
most problematic calculation is determining r generators of a finite index subgroup 
of the Mordell-Weil group, where r is the analytic rank. We used the standard 
Sage method 'E.point_search(height_limit=18,rank_bound=r)', and in case this 
did not come up with enough points we tried some of the remaining curves with 
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'E.gens()'. In several cases these methods did not provide an answer within 48 
hours on a single CPU for a single elliptic curve. For these curves we used the 
method 'MordellWeilShaInformation()' in MAGMA [S], which could handle all our 
problematic curves in a few seconds each. 

The second problematic calculation in the actual code is computing the image 
of coker77Q in K{S,5). We try to factor ideals of K, which are generated by 
elements of possibly very big norm. For example, the curve E^, for d = 1/94, has 
analytic rank 1 and the numerator and denominator of the image of the point of 
infinite order in K{S, 5) have about 600 digits and Sage was not able to factor 
the corresponding ideal. As we already knew that the image is trivial, since the 
dimension of cokerr^q was zero, we could skip this calculation. Considering this 
information in the algorithm all curves we tried worked fine. This problem might be 
avoidable by trying another strategy working modulo primes. The rest of step 1 is 
no problem for moderately chosen d — u/v, since it is mainly prime factorization of 
integers and of rational polynomials of degree 25 (to divide points by 5), as well as 
calculating isogenics and analytic ranks. On a desktop computer one could produce 
in a few hours a database of a few thousand curves. 

Remark 6.1. Step and 2 do not use any assumptions, but for step 1 we assume 
the Birch and Swinnerton-Dyer conjecture in case the elliptic curve is of analytic 
rank r > 2, to conclude that the calculated sublattices A and A' are of finite index 
and that the Tate-Shafarevich groups are finite. Thus, only in case that both elliptic 
curves Ed^ have analytic rank r < 1 the result of the algorithm about ^ni(Bdi^d2) 
modulo squares is completely unconditional. 

7. Results 

Given the above described algorithm, one can produce in short time millions 
of examples of abelian surfaces over Q, such that either the order of the Tate- 
Shafarevich group is a square or five times a square. In case the two elliptic curves 
were both of analytic rank r < 1 these examples are completely unconditional. We 
constructed two databases of elliptic curves using step 1 of the algorithm. The first 
database consists of all elliptic curves Ed, d = u/v, with u,v positive integers and 
max(u,u) < 50,000, such that the conductor of Ed is bounded by C = 10®. The 
second database consists of all elliptic curves, such that max(u, < 100, for u,v 
positive integers. 

Database 1 contains 2212 elliptic curves, all of them having analytic rank r < 
2. It is likely that there are no further elliptic curves of conductor at most 10®, 
which have a rational torsion point of order 5, since there is no such curve with 
4617 < max(u,i;) < 50,000. The database is described in more detail in Table 1, 
where we state for each analytic rank the number of elliptic curves with conductor 
at most 10® and with max(M, v) < N. Database 2 contains 6, 087 elliptic curves. 
All of them have analytic rank r < 3. See Table 2 for more details. In the following 
we will present the results of step 2 of the algorithm applied to the two databases 
described above. 

Database 1 yields 2,445,366 abelian surfaces Bd^^d^- It turns out that 47.01% 
of these surfaces have a Tate-Shafarevich group of non-square order. Database 2 
leads to 18,522,741 abelian surfaces. The percentage of the non-square case is 
49.31. The intersection of the two databases consists of 1,391 curves, hence we 
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considered 966, 745 surfaces twice. In total this gives 20, 001, 362 surfaces, of which 
49.16% have a Tate-Shafarevich group of non-square order. 



N 


#Ea, C = lO'' 


#{r = 0} 


#{r = 1} 


#{r = 2} 


50,000 


2,212 


987 


1,109 


116 


4,617 


2,212 


987 


1,109 


116 


3,375 


2,211 


986 


1,109 


116 


3,072 


2,210 


986 


1,108 


116 


2,695 


2,209 


986 


1,107 


116 


2,000 


2,200 


982 


1,102 


116 


1,000 


2,174 


963 


1,095 


116 


900 


2,170 


961 


1,093 


116 


800 


2,159 


956 


1,088 


115 


700 


2,145 


951 


1,079 


115 


600 


2,119 


941 


1,063 


115 


500 


2,088 


921 


1,052 


115 


400 


2,066 


912 


1,039 


115 


300 


1,993 


872 


1,009 


112 


200 


1,818 


786 


929 


103 


100 


1,391 


616 


697 


78 


50 


845 


394 


405 


46 



Table 1. Database 1: Number of elliptic curves with conduc- 
tor < 10^ and max(u, v) < N. 



N 


#Ed 


#{r = 0} 


#{r = 1} 


#{r = 2} 


#{r = 3} 


100 


6,087 


2,390 


3,038 


633 


26 


90 


4,959 


1,987 


2,463 


490 


19 


80 


3,931 


1,597 


1,940 


380 


14 


70 


2,987 


1,235 


1,455 


287 


10 


60 


2,203 


925 


1,074 


198 


6 


50 


1,547 


660 


760 


123 


4 


40 


979 


412 


494 


70 


3 


30 


555 


245 


277 


33 




20 


255 


130 


115 


10 




10 


63 


40 


22 


1 





Table 2. Database 2: Number of elliptic curves E4, such that 



max(u, v) < N. 



We did two different experiments with the two databases. In experiment 1 we 
investigated how the rank influences the squareness of the Tate-Shafarevich group. 
We Hst the result in Table 3 for Database 1 and in Table 4 for Database 2. The 
first three, respectively four, entries correspond to pairs {Ei,E2) with the same 
analytic rank. The following three, respectively six, lines correspond to pairs with 
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different analytic ranks and the final line corresponds to pairs with analytic rank 
r < 1. If we consider abelian surfaces of fixed analytic rank of at least 4 then the 
density of the surfaces with square Tate-Shafarevich group seems to be significant 
larger than 0.5. However the surfaces with rank larger than 2 inside our family are 
conjectured to have density zero and our database contains very few such cases. 
The calculations with curves of rank r < 1 all show that the non-square case 
happens in about 50% of all cases. For both experiments we list how many abelian 
surfaces -B^j occur in each of the cases, we state the percentage of the surfaces 
with square Tate-Shafarevich group, and we give the percentage of in how many 
cases the parity of the rank of the abelian surface agrees with the parity of the 
exponent of the regulator quotient (RE). Note that the results are unconditional 
in case ik{Ei) < 1, for both If one of the analytic ranks is at least 2 then we 
assume that it is at least the Mordell-Weil rank. 





C = 10«, N = 50,000 






%m □ 


% RE = rk (2) 


r = 
r = 1 
r = 2 


486,591 
614,386 
6,670 


54.041 
58.614 
92.039 


100.00 
63.51 
55.53 


r = 0, r = 1 
r = 0,r = 2 
r = l,r = 2 


1,094,583 
114,492 
128,644 


46.634 
52.867 
74.314 


83.44 
47.96 
42.48 


r < 1 


2,195,560 


51.628 


81.53 



Table 3. Results of experiment 1 for Database 1. 





.V = 100 






%m = □ 


% RE = rk (2) 


r 


= 




2,854,855 


48.598 


100.00 


r 


= 1 




4,613,203 


48.882 


80.91 


r 


= 2 




200,028 


73.031 


44.03 


r 


= 3 




325 


98.154 


51.08 


r 


= 0,r 


= 1 


7,260,820 


51.366 


91.02 


r 


= 0,r 


= 2 


1,512,870 


50.567 


71.36 


r 


= 0,r 


= 3 


62,140 


49.891 


52.73 


r 


= l,r 


= 2 


1,923,054 


52.717 


59.50 


r 


= l,r 


= 3 


78,988 


60.632 


46.23 


r 


= 2,r 


= 3 


16,458 


84.470 


48.23 


r 


< 1 




14,728,878 


50.051 


89.59 



Table 4. Results of experiment 1 for Database 2. 



In experiment 2 we looked for the behaviour of the distribution of square and 
non-square order Tate-Shafarevich groups for increasing conductor, respectively 
height, of the elliptic curves. Hence we filtered Database 1 for different values of 



ON THE DENSITY OF ABELIAN SURFACES WITH #m = 5n^ 



15 



conductor bounds C and Database 2 for different values of height bounds N. For 
low bounds, the non-square case was less likely. When we increase these bounds 
this frequency tends to approximately 50%. The results of experiment 2 is given 
in Table 5 for Database 1 and Table 6 for Database 2. Note that for some of the 
surfaces we assume the weak form of the Birch and Swinnerton-Dyer conjecture 
mentioned above. 



c 






%m = □ 


% RE = rk (2) 


1,000,000 


2,212 


2,445,366 


52.990 


77.84 


800,000 


1,966 


1,931,595 


53.232 


77.16 


600,000 


1,683 


1,415,403 


53.758 


76.06 


400,000 


1,351 


911,925 


54.215 


75.24 


200,000 


924 


426,426 


55.001 


73.91 


100,000 


623 


193,753 


57.074 


74.29 


80,000 


547 


149,331 


57.776 


74.03 


60,000 


470 


110,215 


57.990 


72.75 


40,000 


376 


70,500 


59.306 


73.34 


20,000 


245 


29,890 


61.288 


71.72 


10,000 


152 


11,476 


62.182 


72.59 


5,000 


110 


5,995 


59.783 


71.79 


1,000 


45 


990 


65.556 


76.77 



Table 5. Results of experiment 2 for Database 1. 



N 






%m ^ □ 


% RE = rk (2) 


100 


6,087 


18,522,741 


50.694 


84.14 


90 


4,959 


12,293,361 


50.821 


83.66 


80 


3,931 


7,724,415 


50.941 


83.32 


70 


2,987 


4,459,591 


51.235 


82.51 


60 


2,203 


2,425,503 


51.461 


82.00 


50 


1,547 


1,195,831 


52.211 


80.85 


40 


979 


478,731 


52.764 


79.92 


30 


555 


153,735 


54.157 


77.12 


20 


255 


32,385 


56.384 


77.11 


10 


63 


1,953 


67.179 


74.04 



Table 6. Results of experiment 2 for Database 2. 



The two ways we ordered the elliptic curves, via conductor or via height, are 
natural ways of ordering elliptic curves. It is conjectured that densities obtained 
concerning these orderings agree. In both cases the densities seem to exist and are 
around 0.5. This is in contrast to the results of Poonen and StoU (TUI, who showed 
that for genus 2 curves the density of the non-square Jacobians is about 0.13 and 
this density tends to zero, as the genus goes to infinity. 

We end by giving some heuristics why we expect the density to be 50%. We 
expect that for a random pair {di — ui/vi.d^ — U2/V2) in Q* x Q* the global 
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factor is a square for 50% of the abelian surfaces and that the local factor is a 
square for 50% of them, too. We also expect these distributions to be independent. 
Using the 18, 522, 741 pairs obtained from the second database, we get numerical 
evidence for the independence, as illustrated in the following table. 



N = 100 


global quotient = □ 


global quotient / □ 


local quotient = □ 


26.08% 


25.26% 


local quotient ^ □ 


24.04% 


24.61% 



Recall that the exponent of the local quotient equals — #(Ti U T2) + #(C/i H C/2), 
hence one could prove the expected densities for the local quotient by showing that 
the probability that the set (Ti UT2) has an even number of elements is independent 
of the probability that the set {Ui n U2) has an even number of elements. The 
corresponding numerical result for Database 2 is gathered in the following table. 



N = 100 


#(c/i n t/2) = (2) 


#(f/i n [/2) = 1 (2) 


#(ri u T2) = (2) 


46.71% 


1.80% 


#(ri u T2) = 1 (2) 


49.55% 


1.95% 



The global quotient is harder to control. The exponent of the torsion quotient 
equals 3 on a density 1 subset of the pairs {di,d2), see [H Proposition 4.6]. The 
results of Tables 3-6 suggest that the squareness of the regular quotient and hence 
the squareness of the global quotient is not independent of the parity of the rank. 
If both ranks of the elliptic curves Ei and E2 are equal to 0, hence are even, the 
regulator quotient equals 1, hence is a square. If one elliptic curve is of rank and 
the other is of rank 1, then the regulator quotient is a non-square if and only if 
cokerTyq is trivial modulo torsion, where 77 is the usual isogeny belonging to the 
elliptic curve of rank 1. In Database 2 we have the following situation. For the rank 
1 curves it happens in about 91.2% of the cases that ijq is surjective on the free 
part. In case both ranks are equal to 1, the regulator quotient is a square in about 
80.9% of the cases. For the complete second database we get that the parity of the 
exponent of the regulator quotient agrees with the parity of the rank in 84.14%. If 
we consider only all the elliptic curves of rank < 1, then we have that for abelian 
surfaces i?di,(i2 of even rank the regulator quotient is a square in about 88.2% of 
the cases, and for abelian surfaces Bdi.ds of odd rank the regulator quotient is a 
non-square in about 91.0% of the cases, which together gives 89.6% of agreement. 
We have the following situation for the complete Database 2. 



N = 100 


rk(Brf,.rf,) =0 (2) 


Yk{Bd,.d,) = 1 (2) 


regulator quotient = □ 


42.067% 


7.931% 


regulator quotient =/= □ 


7.927% 


42.075% 



In contrast to the global quotient, the squareness of the local quotient seems to 
be independent of the parity of the rank of the abelian surfaces. Again we give the 
numerical results for Database 2. 



N = 100 


rk(Bd„<iJ = 0(2) 


Yk{Bd„d,) = 1 (2) 


local quotient — □ 


25.670% 


25.675% 


local quotient ^ □ 


24.324% 


24.331% 
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